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Abstract 

We construct a generalised expression for the normal ordering of (a + 
a') m for integral values of m and use the result to study the quantum 
anharmonic oscillator problem in the Heisenberg approach. In particular, 
we derive generalised expressions for energy eigen values and frequency 

2 -2 * 

shifts for the Hamiltonian H = ^- + 2_ + —x m . We also derive a closed 
form first order multi scale perturbation theoretic operator solution of this 
Hamiltonian with a view to generalise some recent results of Bender and 
Bettencourt . 
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1 Introduction 

The simple harmonic oscillator (SHO) model can be used to understand a wide 
variety of physical phenomena ranging from problems in Newtonian mechanics 
to those in quantum field theory. This solvable model describes small oscillation 
of a system about the mean position of equilibrium under the action of a linear 
restoring force. However, for a real physical problem one has to incorporate the 
anharmonicity in the model Hamiltonian. In this case the equation of motion is 
not exactly solvable and, in-fact, the anharmonic oscillators (AHO) provide one 
of the simplest examples of quantum mechanical systems which can not be solved 
without making use of approximation techniques. This was why AHO problem 
got considerable attention of the physicists and mathematicians from the very 
beginning of the subject. It is an outstanding curiosity to note that studies in 
the AHO problems have enriched the subject of large order perturbation theory 
[EB, divergent expansion of quantum mechanics Laplace transformation 

representation of energy eigen values ||-Jl2|, computational physics, Pade and 
Borel summation of perturbation series pa]. 



In the present paper we shall study the generalised quantum anharmonic 
oscillator problem in the Heisenberg representation by making use of two simple 
normal ordering theorems derived by us for the expansion of (a + a 1 *)™. We shall 
see that the merit of the present paper is its simplicity. For example, the present 
method provides, on the one hand, a straight forward mathematical frame work 
to construct expressions for the energy eigen values as well as frequency shifts 
and, on the other hand, generalises some of the recent results on the topic by 
Bender and Bettencourt |lj - ^5[ . 

The Hamiltonian of a generalised anharmonic oscillator having unit mass 
and unit frequency is written as 



H = Hn + ±x r ' 



a , x i \_ r, 

2 ^ 2 m J/ 



(1) 



with 



x=-j=(a! + a). (2) 

Here we have chosen to work in units in which % = 1. Equation of motion 
corresponding to (|l|) is 

x + x + Ax™" 1 = 0. (3) 

Equation (|^) can not be solved exactly for m > 2. However, large number of 
approximation methods are available for solving (||) for particular values of m. 
These include perturbation technique [Q, variation of parameters and Taylor 
series approach |fl7|| . Ordinary perturbation technique leads to the unwanted 
secular terms that grows up rapidly with time. Secular terms are unwanted 
because they are in conflict with the physical requirement that the solution 
be finite. There are some methods which are successfully used to sum up the 
secular terms for all orders. Multiscale perturbation theory (MSPT) is one of 
those methods 

Anharmonic oscillator problem can be approached in two different ways in 
quantum mechanics. One way is the c-number approach or the Schroedinger ap- 
proach and the other way is the operator approach or the Heisenberg approach. 
There is an extensive amount of literature on c-number approach p]- |l3] , ^8| but 
the situation is different for the Heisenberg approach because the complicated 
operator algebra tends to pose serious mathematical problems. Our main ob- 
jective in this work is to provide a theoretical frame work which is free from 
calculational difficulties. The first operator solution of a quantum oscillator was 
given by Bender and Bettencourt Jl^,|l^] using MSPT. They generalised the 
existing theory of MSPT into an operator approach to get the zeroth order so- 
lution involving a quantum operator analog of the classical first order frequency 
shift and interpreted it as an operator mass renormalization that expresses the 
first order shift of all energy levels. Somewhat later Mandal ||l9| approached the 
problem (up to first order) using Taylor series method. All these solutions were 
found for the quartic oscillator only. Application of the Taylor series method to 
sextic and octic oscillator problems are now in order |2C[ |. 
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The Schroedinger approach to AHO problem in quantum mechanics has 
made a significant contribution for the development of the perturbative tech- 
niques as we use today but the suitability of the other approach (Heisenberg 
approach) has not yet been studied in detail. Till now we don't have any opera- 
tor solutions of higher anharmonic oscillators although these might be relevant 
to various physical problems. Keeping these in mind, we will suggest a novel 
method in which MSPT results for anharmonic oscillators in general can be 
found. To that end we prove in section 2 two theorems to construct a normal 
ordered expansion of (a + a' ) m . We obtain a generalised expression for the 
energy eigen values in section 3. We devote section 4 to construct a generalised 
solution for the equation of motion (||) and specialise our result to reproduce 
some of the existing results as useful checks on the generalised solution obtained 
by us. 



2 Operator Ordering theorems 

On a very general ground one knows that in quantum mechanics and quantum 
field theory proper ordering of the operators plays a crucial role. Let /(a, a^) 
be an arbitrary operator function of the usual bosonic annihilation and creation 
operators a and a* which satisfy the commutation relation 

[a,a f ] = l. (4) 

One can write f(a, a') in such a way that all powers of always appear to the 
left of all powers of a. Then f(a, a') is said to be normal ordered. In this paper 
we want to write (a + a t ) m in the normal ordered form for integral values of 
m. Traditionally for a given value of m this is achieved by using a very lengthy 
procedure which involves repeated application of (^). One of our objectives in 
this work is to construct a normal ordered expansion of (a + a^) m without taking 
recourse to such repeated applications. 

We shall denote normal ordered form of / by /at- On the other hand : / : 
will denote an operator obtained from / by arranging all powers of a) to the 
left of all powers of a without making use of the commutation relation in (f|) . 
Now if / = aa^ then = a^a + 1 and : / := a^a. Therefore, we can write 

: (aW)™:— (a^ + a) m := a m + m C 1 a^a m - 1 + ...+ m C r a r 'a m - r + .. + a rn . (5) 

Thus in this notation : (a + a^)" 1 : is simply a binomial expansion in which 
powers of the a* are always kept to the left of the powers of the a. To write 
(at + a)jy we shall proceed by using the following theorems. 

2.1 Theorem 1. 

: (at + a) m : (a f + a) =: (a f + a) m+1 : +m : (a f + a)™" 1 : . (6) 
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Proof: From (|) : (at + a) m : (at + a) can be written in the form 

: (at + a) m : (at + a) 

= at m+1 + ( m d + m C )at m a+ ( m C 2 + m d)^ 1 a 2 + .... + ( m C r + m C r -i)a r 

+ ( m Ciat m_1 + 2 m C 2 at m_2 a+ ... + r m C r at"*~V- 1 + ...) . 

_ (7) 
The theorem in (|6j) can be obtained by simply using the following identities 

in the above 



(a r at)7v = a^a r + ra r l , 

r n C r = n n - l C r -i, 

and m+1 C r+ i = ( m CV+ m CV+i). 



(8) 



Note that first identity in (||) can be proved with the help of the general operator 
ordering theorems while the other two are trivial. 



2.2 Theorem 2 : 

For any integral values of m 

(at 



with 



r=0,2,4.. 

(r-1)! 



and 



2 (S-i)(|_i)! 
to = t 2 = 1 



l C r : (at + a) r 



/or r > 4 



Proof: 

Using theorem 1, we can write 



(at + a) N 
(at + a)% 
(at + a)% 
(at+a) 4 w 
(at + a)^ 



(at 
(at 
(at 

(at + a) 

(at + a) 5 : +10 : (at + a) 3 : +15 : (a f + a) 



(at + a) : 

(at + a) 2 : +3 



(9) 

(10) 
(11) 



(12) 



(at + a )% = : (at + a) 9 : +36 : (at + a) 7 : +378 : (at + a) 5 
+ 1260 : (at + a) 3 : +945 : (at + a) : . 



From ( |12| ) we venture to identify the general form of the above expansion for a 
given value of m as 

m 

(o t + o)^= t r m C r : (a? + a) m - r : . (13) 

r=0,2,4.. 
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It is easy to check that (|l3|) gives all the expansions of ( |T^ ) such that ( |l~3| ) is 
true for m — 1,2,.., 9. The general validity of ( p~3[) can be ensured by using 
the method of induction. This establishes that (|13[) gives the normal order 
expansion of (at + a) m for any arbitrary integer to. We shall now use this 
normal ordered expansion to study the anharmonic oscillator problem. 



3 Energy eigen values 

The first order energy eigen value Ei is < n\H\n >, where \n > is the number 
state. In the literature there are two lengthy procedures [^| to obtain E\. The 
first one is the usual normal ordering method. This method involves iterative 
use of @j and number of iteration increases very first as to increases. For a 
given value of to we need [2 m — (to + 1)] iterations. This shows that for large to 
the construction of the expansion (at + a)^ becomes formidable. In the second 
procedure one proceeds by making repeated application of x operator on the 
number state. As in the normal ordering method this procedure is also equally 
lengthy. In the following we implement theorem 2 to derive an uncomplicated 
method to construct an expansion for E\ for the Hamiltonian (|]). 

Using the result in @ we can write the expression for E\ for any integer 
value of to in the form 

Ei = (n + \) + ^- < n\ ££=0,2,4.. tr m C r : (at + a) m - : |n > 

- (n+|) + ^- < n|E™ 2 4 U m C r m - r C m - r a^a^\n> 

/„ i 1 \ i A -f ms~i m — r/~< n/^i I rn — r \ \ 

- ( n + 2 > + 2^^=0.2.4..^ Um-r )'■ 

" " (14) 

The expression (14) involves summations which are easy to evaluate and thereby 
avoids the difficulties associated with earlier iterative procedures. Although 
somewhat forced, it may be tempting to compare the simplicity sought in our 
approach with the use of logarithm in a numerical calculation or use of integral 
transforms in solving a partial differential equation. It is of interest to note 
that expression similar to ( fil} ) can also be constructed for the first order energy 
eigen value of a Hamiltonian in which the anharmonic term is a polynomial in 
x. This type of Hamiltonians are very frequently used in nonlinear optics. 



4 MSPT solution of the generalised quantum an- 
harmonic oscillator. 

We want to obtain MSPT operator solution of equation (||) for arbitrary inte- 
gral values of to. The essential idea behind our approach is that the quantum 
operator analog of the classical first order frequency shift is an operator function 
( fl(Ho) ) of the unperturbed Hamiltonian Hq and secondly a correct solution 
should reproduce the first order energy spectrum. Now the general form of the 
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zeroth order solution of (^) is given by 

x (t) 



^ [x(0)cos(t + \fl(H )t) + cos(t + Xfl(H o )t)x(0) 
x(0)sin(t + Xn(Ho)t) + sin{t + XQ(H o )t)x(0)} 



(15) 



where G(n) is a normalization factor. So our task is to find out fl(H ) and 
G(n) in general. From equation ( |l4|) we obtain the energy difference for two 
consecutive energy levels as 



Un,n-1 — {Ei) m ,n ~ (-El)m,n— 1 

= 1 + Xu>(m, n) 



1 



x 



E(m-2) , m( ~i 
r=0,2,4.. b r 



2"3" m 



mr~i m — r 



c, 



(16) 



Since the correct quantum operator solution has to give the first order energy 
spectrum, we should have 

< n—l\xo(t)\n >=< n— l|xo(0)|n > cos[t+Xu)(m, n)t]+ < n—l\po(Q)\n > sin[t+Xuj(m, n)t] 

_ _ (17) 

Equation ( |l5|) and ( |l7|) impose restrictions on our unknown functions SI(Hq) 
and G(n). The condition imposed on fl(Ho) is 



or. 



< n\Q(H )\n > + < n - l\Q(H )\n - 1 >= 2w(m,n) 



Q(n + ~) + fi(n - ~) = 2w(ro, n) 



(18) 



(19) 



For a particular m right hand side is a known polynomial in n and our job is 
simply to find out fi(n + |). We obtain this as 



n(n+-)=2 



^(-l)"" fc o;(m,fe) 



+ (-1) 



t r . 



7T1-2 

2 2 m 



(20) 



Substituting the functional form of 0(n + |) or Q(Hq) in ( |l5| ) if we impose 
condition (llTn we will get 



G(n) = 2cos 



^(n(„ + i)-n(„_r 



(21) 



The results in ( |2(i| ) and ( |2l] ) when substituted in (JlJ) solve the generalised 
quantum anharmonic oscillator problem. 



4.1 Specific results and their comparison with the existing 
spectra: 

Although we have solved the quantum anharmonic oscillator in general, it is not 
possible to compare our solution directly with other existing results since the 
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present study happens to be the first operator solution of the generalised anhar- 
monic oscillator. There are methods to calculate first order classical frequency 
shift for particular m in the appropriate classical limit (x(0) = a, i(0) = 0). 
Here we calculate some specific results from our general expressions and com- 
pare them with the existing results. 
For m = 4 we have, 

_ . 1 . 3n 3 3 . 1, . . 

n{n+ 2 ) = ~ + 8 = 4 (n+ 2 } ' (22) 



Therefore, 

M{Ho) = 4-ffo ,. )0 . 
and G(n) = 2cos(^) . ( A> 

In terms of (p3[) the total solution for the quantum quartic anharmonic oscillator 



*(*)L=4 = 2cos (M*) [x(0)cos[t + ^H ] + cos[t + ^H o }x(0) 
+ x(0)sin[t + ^fH o ] + sin[t+^fH o ]x(0)] 



This exactly coincides with the solutions given in ||14|,|15j and |19| and gives the 
correct classical frequency shift J2j| in the limit x(0) = a, i(0) = 0. Similarly 
we have 

*(*)U=6 - 2cos( k n) [x(0)cos[t + f {Hi + \)t] + cos[t + f (Hi + j)t]g(0) 
+ i(0)sm[i+ f (ff n 2 + i)t] + sm[i+f (Hi + \)t]x(0)] 

(25) 



^(t)U=8 = 2co8[ ^; 6n3+3)t] |>(0)aw[f + f^o 3 + 5H )t] + cos[t + f^(4ff 3 + 5H )t}x(0) 
+ x(0)sin[t + ^(4i?3 + 5Ho ) t ] + sin [ t + ^(4i/3 + 5H )t]x(0)] 

(26) 

and 

*(*)U=io = 2coa[ 6 ¥ i 1 (n 3 +2ra)] Now* + fH# 4 + i g o + 

+ «»[* + ^(i/ 4 + |ff 2 + £)f]s(0) + ±(0)*in[t + ^(i/ 4 + \Hl + ±)t] 
+ sln lt+l§^(H£ + lHZ + ^)t}±(0)] ■ 

(27) 

The solutions for sextic and octic oscillator exactly coincides with the solution 
obtained by us using Taylor series approach |2^] and all the solutions give the 
correct classical frequency shifts in the appropriate limit [ |23J2^| . 



5 Summary and concluding remarks 

We conclude by noting that depending on the nature of nonlinearity in a physi- 
cal problem the treatment of higher anharmonic oscillators assumes significance. 
But studies in such oscillators (for m > 4) are not undertaken in the Heisenberg 
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approach presumably because the existing methods tends to introduce inordi- 
nate mathematical complications in a detailed study. In the present work we 
contemplate to circumvent them by proving a theorem for the expansion of 
(at + a) 7 }}- Thus the results of the present work are expected to serve a use- 
ful purpose for physicists working in nonlinear mechanics, molecular physics, 
quantum optics and quantum field theory. 
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